A SURVEY OF DYNAMICAL MATRICES THEORY 

LIN ZHANG AND JUNDE WU 

Abstract. In this note, we survey some elementary theorems and proofs concerning dynamical matrices 
theory. Some mathematical concepts and results involved in quantum information theory are reviewed. 
A little new result on the matrix representation of quantum operation are obtained. And best separable 
approximation for quantum operations is presented. 
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1 . Introduction and preliminaries 

Positive linear maps on some operator algebras are a very important subject of both the mathemat- 
ical and the physical literature for several years. The images of positive operators acting on a given 
Hilbert space under such a map are positive operators acting on the same Hilbert space. A map <1> is 
called ^-positive for some k e N if the tensor product CD ® Id* is positive. We call <& is a completely 
positive (CP) when it is ^-positive for any k e N. Completely positive maps (CP maps, for short) de- 
scribe the dynamics of open quantum systems. The structure of the set of CP maps is well understood 
due to the theorems of Stinespring [12], Kraus [8], and Choi O. Choi's theorem is also proved by 
another simple approach in [11]. 



In this paper, only finite dimensional complex vector spaces are considered. An column vector in a 
complex vector space is denoted by \<f>), the symbol <p is a label, while |-) denotes that the object is a 
complex column vector. This notation for complex vectors is called Dirac notation. Throughout the 
paper, t,t and * stand for Hermitian conjugate, transposition and complex conjugate, respectively, of 
matrices with respect to a given orthonormal basis. Given a vector 

10) - [0i, 
its dual is defined as 

(<f>\ = [<p\ <f>l ■■■ (/>*] = (|0»t. 
Given the vectors \<p), \<p), the inner product between two vectors is denoted by (4>\<p), which is defined 
as follows: 

The norm of a vector \(p) is defined as ||0|| = ^((p\(p). Unite vectors are those vectors with unit norm. 
Two vectors are orthogonal it they have zero product. The outer product of the given vectors \<p) and 
is given by 



01 

4>2 

<Pd 



W\ <P*2 ' ' ' V*d\ = 



<f>np\ 4>i<f* 2 

<!>2<P\ 4>2<P*2 
</>d<P*i <f>d<P* ? 



tWd 

0d<P*, 



A set of vectors {|vyt>}^ =1 in a vector space *V is orthonormal if the vectors are normalized and orthog- 
onal, that is, {vj\vj) = Sij. If, in addition, n = dim'V, this set of vectors form an orthonormal basis 
for < V. Here we have a simple but useful fact that Y!k=i \ v k)(vk\ = In for given an orthonormal basis 
{|vfc)}^ =1 in a vector space < V. This called the completeness relation. 



Quantum states will now be introduced. A quantum system is a physical system that obeys the 
laws of quantum mechanics. Let us assume that we are given two quantum systems. The first one is 
owned by Alice, and the second one by Bob. The physical states of Alice's system may be described 
by states in a Hilbert space 'Ha of dimension cIa = N, and in Bob's system in a Hilbert space 'Hb of 
dimension ds = M. The tensor product is a ubiquitous mathematical operation which can be used 
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to combine vector spaces to form a larger vector space. Given two vector spaces *V and <W ', we can 
combine them to form the vector space r V®'W, with dim('V® < H / ) - dim('V)xdim('ll / ). The bipartite 
quantum system is then described by vectors in the tensor-product of the two spaces - Ha ® H B , 
and dimCK) = cIacIb- A pure state of dimension d can be represented by a J-dimensional complex 
unit vector \if/). For real 9, the vectors \if/) and e' e \ft) represent the same state. More generally, a d- 
dimensional quantum state is represented by a d x d complex matrix p, also called a density matrix, 
which is a non-negative linear operator, acting on a complex Hilbert space fi, with trace 1. A pure 
state can be represented either by its state vector \if/), or by its density matrix p = States which 

are not pure are called mixed states. A simple test for whether a state p is pure or mixed is to take the 
trace of p 2 : tr(p 2 ) = 1 if p is pure and tr(p 2 ) < 1 if p is mixed. A mixed state can be expressed as a 
mixture of pure states in many different ways. 



Suppose that |v> € *V, |w) e < W. The vector |v) ® \w) e *V ® < W. The vector |v> ® \w) is computed as 
follows: 





wi\v) 




W\ 




Vl 


v) ® \w) - 


w k \v) 


if \w) = 


w k 


and |v) - 


Vk 




. w„\v) 




. w d B . 




. v d A . 



Similarly, the tensor product of two given matrices will be explained as follows: with the orthonormal 
bases {|m)}(m = 1, . . . ,g?a) of "Ha and ||p)}(p - 1, . . . ,dg) of "Hb, respectively, the orthonormal basis 
of can be described as \\m) ® \p) = \mp)}(m - l,...,dA',n - !>•••> ds) (throughout the present 
paper, Roman indices correspond to the subsystem A and Greek indices to the subsystem B.) for which 
two types of ordering are very important such as: 

(i) Ordering of type-I: 

(111), |21), . . . , |d A l>; . . . ; \lfi), |2ju>, . . . , \d A n); . . . ; \ld B ), \2d B ), . . . , \d A d B )\. 

(ii) Ordering of type-II: 

(111), |12), . . . , \ld B ); . . . ; |ml>, |m2), . . . , \md B ); . . . ; \d A l), \d A 2), . . . , \d A d B )}. 

SB^Ha) and &CH B ) means that the set of all bounded linear operators on "H^a and <H B , 
respectively. If I e 3S{ < H A ) and Y € 3§{ < H B ), then X ® F e ^CH). Suppose that the matrix- 
representations X = [x mn ] and 7 = [y^y] for X and F with respect to the given orthonormal bases 
{N)}^ =1 and (|ju)}* j are given, respectively. Then there are several different matrix-representations of 
X ® F with respect to the corresponding orthonormal bases of different orderings. For the ordering of 
type-I, the matrix representation of X ® F is 
y u X y n X ••• v Wb X 



X®Y = 



y 2 \X y 22 X 



y2d B X 



yd B iX yd B iX ■ ■ ■ yd B d B X 
while for the ordering of type-II, the matrix representation of X ® F is 
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X®Y = 



x 2 \Y X22Y 



xu A Y 

X2d A Y 



Xd A \Y x dAl Y ■■■ x dAdA Y 
The ordering of type-I will be employed throughout the present paper if unspecified. For tensor prod- 
uct, we have the following rules: given two matrices S and T acting on vector spaces <V and < W, 
respectively, vectors \x) e *V and |y) e 1V, then 

(S ® T)(\v) ® |w» = (5|x»®(r|y»,Tr(5 ® T) - tr(5)tr(T), (5 ® T) 1 " - S 1 " ® T 1 ". 
If X, Y act also on *V, *W respectively, we have (5 ® ® 7) = ® 7T. Obviously, tensor product 
is a bi-linear map. 



The description of subsystems of a composite quantum system is provided by the reduced density 
operator, which is so useful as to be virtually indispensable in the analysis of composite quantum 
systems. Suppose we have physical systems A and B, whose state is described by a density operator 
Pab- The state space of the composite quantum system AB is denoted by S'i'H), similarly, SH^Ha) for 
subsystem A and Q}^Hb) for subsystem B. The reduced density operator for system A is defined by 
^b(Pab) = Pa, where Tr# is a map of operators know as the partial trace over system B. The partial 
trace is defined by 

Tr fl (|a 1 X«2l ® \bi)(b 2 \) = |ai><fl2|tr(|fei><62l), 
where \a\) and |a2> are any two vectors in the state space of A, and \b\) and I&2) are any two vectors in 
the state space of B. The trace operation appearing on the right hand side is the usual trace operation 
for system B, so tr(|£iX&2l) - (^2^1)- In fact, Tr A = tr®Idg, Trg = Id^^tr and Tr = tr®tr. We have de- 
fined the partial trace operation only on a special subclass of operators on AB. More generally, for any 
matrix Z acting on *Ha®7~(b, we have a block construction on Z: Z = [Z^ijx, v - \,...,ds = dim "Kg), 
where each Z MV is a scalar matrix of size dA x g?a(^a = dim "Ha)- Therefore 

Z = ^ =1 ^ v ®|//Xv|. 
Now the partial trace over system A is provided by 

Tr A (Z) = Zj v=1 tr(Z MV )\p)(v\ = [tr(Z^ v )], 
while the partial trace over system B is given by 

Tr B (Z) = Ej^j Z^trdjuXvl) = I % Z m . 
The partial trace over the composite quantum system AB is 

Tr(Z) = 2jj tr(Z w ). 

The quantum operations formalism is a general tool for describing the evolution of quantum sys- 
tems in a wide variety of circumstances, including stochastic changes to quantum states. A simple 
example of a state change in quantum mechanics is the unitary evolution experienced by a closed 
quantum system. The final state of the system is related to the initial state by a unitary transformation 

U, 

p _» £(p) = UpUK 
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Unitary evolution is not the most general type of state change possible in quantum mechanics. Other 
state changes, described without unitary transformations, arise when a quantum system is coupled to 
an environment or when a measurement is performed on the system. This formalism is described in 
detail by Kraus. In this formalism there is an input state and an output state, which are connected by 
a map 

£(p) 
p ~^ tr[£( P )] ■ 

This map is determined by a quantum operation £, a linear, trace-decreasing map that preserves posi- 
tivity. The trace in the denominator is included in order to preserve the trace condition tr(p) = 1. The 
most general form for £ that is physically reasonable, can be shown to be 

£(p) = Z,r>r}. 

The system operators Tj , which must satisfy Yij r^rt < /, completely specify the quantum opera- 
tion. Formally, every quantum operation has to be described mathematically by a completely positive 
complex-linear mapping £, which satisfies tr(£(p)) < 1 for all state p. A quantum operation is called 
quantum channel if it is trace-preserving. 

Given quantum operation £, £4, and £g on corresponding bipartite quantum system with subsys- 
tems A and B, subsystems A, and B, respectively, owing to Jamiolkowski isomorphism, the notion 
of entanglement can be extended from quantum states to quantum operations. A quantum operation 
acting on two subsystems is said to be separable if its action can be expressed in the Kraus form 

£(0 = £*(a£®a*)-(A>a*)1\ 

where and A^ are operators acting on each subsystem and they satisfy that 

E*(A£ ® Af)t(A£ ® A*) < I A ® I j}. 
Otherwise, it is entangled. When the equality is valid, there is a concept of separable quantum chan- 
nel. 



2. Vectorization and realignment of matrices 



Definition 2.1. Representation of matrices as vectors on a higher dimensional Hilbert space is called 
vectorization. It transforms a pXq matrix G into pq x 1 column vector denoted by |G)\ this is done 
by ordering matrix elements, i.e., by stacking the columns of G to form a vector: for example, with a 
px q matrix G = [ga], \G)) is described as 



|G» = 



G(-,l) 



G(;q) 



, where G(-,j) = 



8\j 



U = 1, ■•-,<?)• 



That is, G(-,j) is the jth column vector of matrix G. Dually, ((G| is a 1 x pq row vector defined as 
(|G»)+, i.e., «G| - (|G»)+. (see (1) 



Remark 2.2. (i) Vectorization is obviously linear: for matrices S k and scalars At, 

\Yjk^kSk)) - Yjk^k\Sk))- 
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(ii) Vectorization is inner-product-preserving; i.e. isometry. The Hilbert-Schmidt inner product is 
equivalent to the usual Euclidean inner product of vectors: for square matrices S, T of the same 
size, (S,T) - tr(S^T) - ((S\T)). It is easily shown that vectorization is one-one and onto. 
Therefore vectorization is a unitary transformation from Hilbert-Schmidt matrix space to Hilbert 
vector space. 

(iii) Vectorization is intrinsically related to the tensor product. Consider a square matrix of size 
px p, representing an operator acting on the /^-dimensional Hilbert space 'K. Let {|j')}y =1 be the 
orthonormal basis of 'K for which \j) is column vector with all entries except for jth entry 1. A 
matrix T = [tij] = 2f - =1 UjEij, where E,j = is transformed to the vector 

i,y=l 1,7=1 1,7=1 7=1 (=1 

P P P 

7=1 i=l 7=1 

P 

(2.1) = (T ® I p )\ J] E n )) -{T ® I p )\I p )) = (I p ® T^lp)). 

7=1 

Thus it follows from the identity above that, for any matrices Q,X and R of the same size px p, 
\QXR)) = (QXR) ® I p \I p )) = (Q® I p )(X ® /,)[(/? ® / p )|/ p »] 

= (G ® /„)(* ® / p )[(/ p ® /?*)|/p»] = ((2 ® /„)[(* ® / p )(/ p ® fl 1 )]^)) 
= (G ® W„ ® ® i P Wp)) = (G ® W„ ® 

(2.2) - G®^l^» 
and 

(2.3) |XT» = (X ® 7 p )|r» - (I p ® F 1 )^)). 

(iv) For any matrix F, 

(2.4) «n = on))* = (|y>»*t - (in)/. 

(v) For 5 € 3§(^Ha) and T e ^CHb), where "Ka - *Kb are <i-dimensional Hilbert spaces. For the 
matrix representations S = [sjj] and T = [tij](i,j - 1, . . . , d), we have trg(|S))((r|) = ST^ and 
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tr A (|S»«r|) = S t T*. Indeed, 

d d 

tr B (|5»«r|) - s mn tl v \r B {\mn){nv\) = J] W^tr fl (|m><^| ® |n><v|) 

m,n,/i,v=l m,n,ii,v=\ 
d d 
Y s mnt* MV Snv\m)(fl\ = ^ S mn t* n \m)(n\ 
m,n,/i,v=l m,n,/i=l 
d d d d 

n=\ m=l ju=l n=l 



- 2 1 S|i»><n|Tt =lS r t . 



n=l 

The other identity goes similarly. 

Definition 2.3. Let Z be an <ig x <ig block matrix with each entry of size cIa x <s?a; i.e. Z = [Z^ v ] 
represent an operator acting on Ha ® ^Hg- We define a realigned matrix ^(Z), acting from *Hg ® T/g 
to Ha ® *Wa> °f s i ze ^ x tnat contains the same elements as Z but in different position as 

K(Z) - [|Z„», . . . , |Z dsl »; . . . ; |Z Wa », . . . , IZ^))]. 
In fact, ??(Z) Similarly, we can also define another alignment as ??'(Z) Note 

/iv nv jfv //m 

that alignment of matrices is a one-one linear mapping from the matrix space Md A d B xd A d B (C) onto the 
matrix space M^ X( fi(C). 

Proposition 2.4. For a tensor matrix X ® F w/tfi factor matrix X of size cIa X <5?a <z«^ ^ factor 
matrix Y = \y MV ] of size x ds, Z = [y^X] = [Z MV \. We have: 

(2.5) K(x®y) = |x»«n 

Moreover, a nonzero matrix Z can befactorized as X®Y if and only if rank['R(Z)] = 1. 
Proof. 

K(X®Y) = [\ynX)),...,\yd B iX)y,...;\yu B X)),...,\yd B d B X))] 
= [yn\X)), . . . ,y dB i\X)); . . . ;y ldB \X)), . . . ,y dsdB \X))] 

= \x))\y lu . ..,y dB v,... ;yu B , . . . ,y dBdB \ = WW))) 1 
= ix»«n. 



□ 



For a general block matrix Z, it holds that 

d B d B 

<R(Z) = K{ 2 z ^v ® lju)<v|) = 2 ft&vy ® U*Xv|) 
(2.6) = J] IZ^XI/zv)) 1 - 2 l Z /">XM 
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Before the properties of realignment derived, we need to know one useful operator called swap oper- 
ator, denned as S = J^f acting on "Hn ® Then by simple computations, we have: 

I, J 1 

Proposition 2.5. For any X and Y of the same size N X N. We have: 

( i) S is self-adjoint, unitary, symmetric, and orthogonal; 

(ii) \X t )) = S\X)),L T = S; 
(Hi) S(X®Y)S =Y®X. 

Definition 2.6. With 5 as above, the flip transformation of matrices over a bipartite quantum system 
is defined as 

F(Z) = SZS with T(Z) m , = Z m . 



nv vn 



Similarly, we can define two partial flips as T r (Z) - S Z with T r (Z)>w = Z^n and T C (Z) - ZS with T~ c (Z)«v 
Ziw (where 'r' and 'c' mean that row and column, respectively). Later, we will see that Lj- - S ®S . 

Lemma 2.7. Given any two square matrices X, Y of the same size, we have the following equa- 
tion: 

(2.7) |x®r» = (/®s®/)|x»|y». 

Proposition 2.8. (i) IfX, Y are matrices of the same size N x N, then 

(2.8) |ft(X®y)» = |X»|y»; 
i.e., the vectorization of the matrix \X))((Y*\ is \X))\Y)). 

(ii) Let Zbea matrix of size N 2 x N 2 . Then: \R(Z))) = I ® 5 ® I\Z)), thus Lr = I ® 5 ® /. 

(Hi) lf£l{-) = ® lOO'l) • ® I), then: for any matrices X, Y of the same size N x N, 

I, J I 

(2.9) Q(|X»«Y|) = X® Y* and Q.(X ® Y*) = \X))((Y\ = K(X ® Y*). 

More generally, we have Q(Z) = ( R(Z)for any matrix Z of size N 2 X N 2 . 

Proof, (i) and (ii) follow easily from Lemma l277l 
(iii) Together with Lemma 12771 it follows from (i) that 

|H(x® r)» 

N N 

= I £ (/ ® \i)(j\)X ® Y*(\i)(j\ ® /)» = £ |(7 ® \i)(j\)X ® Y*(\i)(j\ ® /)» 

i,j=\ i,j=l 
N N 

= J (/ ® I0O1) ® CIjXiI ® D\x ® y* » = £ (/ ® \ij)(ji\ ® i)\x ® r » 

1,7=1 i,7=l 

= (/ ® s ® /)ix ® r » = |x»|f*» = ® r )». 
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Hence Q(X®Y*) = ft(X®Y*) = By simple computations, we have also fi(|X»«F|) = X®Y*. 

Since (m/u\ft(Z)\nv) - ft(Z)m^ = Zmn and 

Zy =i (/ ® lOO'DZdOO'l ® /)l»v> - Zy=i S^njimmiv) = (mn\Z\fiv) = Zmn, 
i.e., Q(Z) = ft(Z). In such a way, we obtain the explicit expression for the realignment transformation: 



(2.10) 



K(Z)= 2(/®|i)01)Z(|i)<7l®/) 

(J=l 



for any matrix Z of size N 2 x N 2 . □ 

Next the relationship among f/ie realignment, the transposition, and the flip over a bipartite quantum 
system will be discussed. First recall that the transposition T over bipartite quantum system Pi a ® "Hs 
are defined as 7~(Z) = T"a ® Tg(Z) with T{Z)m^ - Z^v , where 7"a and 7"b are the transpositions with 
respect to subsystems A and B, respectively. Apparently, Ta(Z) and Tb(Z) m// — Zmv . 

nv my nv njU 

Proposition 2.9. ( jj 7", ft and f all are involution; i.e., 7~ o 7~ = ft o ft = f o f = Id. 

(ii) f of - T o f, T o ft ^of and T o ft ± ft o f, where o stands for the composite of 

transformations. 
(Hi) T o<R = <RoT andftoT = T oft. 

(iv) ft' = ToftoT = ToftoT. 

(v) T r = ftoT A oft and r c = ftoT B oft. 



Proof. It is trivially by some computations. For example, [T ° ft(X)] n 



Xiiv and 



[ft o T(X)]n,K = [ft(X)]„m = X?, ; i.e., [T o ft(X)],*, = [ft o 7"(X)]^ which means that T oft-ftof. 

nv vn mn nv nv 

□ 



Others go similarly. 



3. Dynamical matrices for quantum operations 



A density matrix 
(3.1) 



P = 



Pn 



Pd B \ 



of size ds x <i# may be treated as a vector 

p(; 1) 



(3.2) 



lp» = 



p(-,d B ) 



Pld B 



Pd B d B 



= [P»v] 



, where p(-,v) - 



Pi> 



Pd B V 



(V = l,...,fi? B ). 



Suppose that p and cr act on 'Pis and "Ha, respectively. The action of a linear super-operator O : p 

a x d 2 B '- 



cr = (I>(p) = [cr m „] may thus be represented by a matrix L<d = L of size d 2 .xd 2 - 



(3.3) 



|cr» = |0(p)» = L|p» or cr m „ - ^ Lpj^. 
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It can be written concretely as the equation of multiplicity of a supermatrix and a supervector. 





" <r(-,l) 




L n ■ 


L\ v 


(3.4) 


<x(-,ra) 




L n \ 


L n v 




cr(-,d A ) 






Ld A v 



where 



(3.5) 



L m - [Lmn\ 



Ll n 

lv 



Ld A n 

lv 



L ld R 



-inds 



Ld A d B 



L l n 

dav 



Ld A n 

d B v 



p(-,l) 

P(;v) 
p(-,d B ) 



One must be caution here that n and v stand for the block row index and the block column index, 
respectively; while m and p stand for row index and column index of each block. Now we give a 
simple example for a qubit map for later use as follows: 



(3-6) 



en 

C"21 
0-12 
0"22 



Ln 
11 


Ln 
21 


Ln 
12 


Ln 

22 


L 2l 
n 


L21 
21 


L21 
12 


L21 
22 


L 12 

n 


L\2 
21 


L\2 
12 


L\2 

22 


L22 
11 


L22 
21 


L22 
12 


L22 
22 J 



Pll 
P21 
P12 
P22 



Theorem 3.1. Tfte requirement that the image <r is a density matrix, so it is Hermitian, positive with 
unit trace, impose constraints on the matrix L: 

(i) Cr^=0- => L*„ - Lnm. 

(ii) a > =^ [Z^ v= i ^jPpv] > for any state p = [p^]. 
(li/J ?r(o") - 1 => Z^ =1 £™ - <V 

Proof, (i) Step 1: For state p = \y)(y\(y e {1, . . . , d B }), 

(3.7) p^y = <p|p|v> = (p\y)(y\v) = d^Syy. 

Then 



(3-8) 



cr„ 



11 a 



LnmSuyS 



Linn . 

yy 



Since o" 1 " = 0-, it implies that cr mn = cr nm ; i.e., L™ = L* m (r e {1, . . . , d B }). 
Step 2: Setting 



(3.9) 
we have 
(3.10) 



P = |[|arXar| + WW + l<*></8| + |0><ar|](ar,j8 = h . . . ,d B ; a * fi), 
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Hence 

(3.11) @~mn ~ ~z\_L" m ~t~ Lmn + Lnm + Z^mn], 

2 aa PP af} jia 

From the equation cr mn - cr* im , we know 

Lmn + Lmn + Lmn + Z^mo — L nm H~ L nm + Z/ wm + ^raw s 

>3j8 or/3 /3a ora /3jg a-/3 /3a 

i.e., 



(3.12) Z/mn + Lmn — L nm + L nm . 

a/3 pa a p p a 

Step 3: Letting 

(3.13) p = j[\a)(a\ + + V^l|a></3| - V^ipXal], 



we have 



(3-14) P^y - -[(^(Att + J^vyS + V-Uf^Syp ~ ^ ~\8 



Hence 
(3.15) 

which implies that 
(3.16) 

This gives rise to: 



(J~ mtl — — \Lmn + Lmn + y — 1 Lmn — 'V — \ Lmn 1 , 

2 L aa pp a p p a 



®~ tlffl r\ \J-' nm L nm y "\L nm H~ *\/ lX/uiuJa 

2 aa pp a p p a 



(3.17) Lmn — Lmn — — L nm + L nm 

aP Pa a p p a 

Combing (I3TT21 with (13~T7T ) gives that L 

^ /3a 

(ii) is trivial. 

(iii) Because tr( cr) = 1, that is, 

(*A ^A rffl 

(3-18) 1 - ^ CTmrn- ^ ^ Lmmppy. 

m=l m=l/i,v=l 

Step 1: Given p = |y)(y|(y e {1, . . . ,<i#}). So p^y = (/u|p|v) = d My S V y. From the equation (13.181 ). we 
have that 

cfo rf B d A 

(3.19) 1 = ^ ^ L^ySyy = ^ L™»(y € |l,...,<f B }). 

m= 1 /J,v= 1 m= 1 

Step 2: From the equation dl!9l , (13.101 ) and (13.181 ). we have that 

4* d A d A d A 

ZT <mm + / Lmm + y Lmm -\- / Lmm 
aa Z—i Bp Z—iaP Z-J /3a 



'4 



PP /—k ap i—l pa 

m=\ m=\ m=\ m=\ 



i.e., 

<4i «U 

(3.20) Y Lmm+Y Lnm = 0. 

l~l ap I—i p a 

777= 1 m= 1 
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Step 3: It follows from the equation (13.131) and (13.141) that 

" d A d A d A d A 

ZJ jtnm -\- ^ Lmm + *\/ — 1 / Lmm — *\/ — 1 / Lmm 
aa Z.J SB *—laB t—l Ba 



i = i 

2 



0B ^—l aB 

m=l m=\ m=l m=l 



i.e., 

d A d A 

(3.21) V L™ - V L™, = 0. 

m=\ m=\ 

From the equations (13.201) and (13.211) . we get 2„ = i - 0(a ^ /?)■ In summary, 2„ = i ^™ = <W □ 

Note that the property (i) of the proposition 3.1. is not the condition of Hermicity, and in general 
the matrix L representing the super-operator <I> is not Hermitian. However, by the definition of matrix 
realignment we can define the dynamical matrix or Choi matrix (see [ 131 [151 ): 

Da> - "RiL) With Dm? = Lmn. 

nv fiv 

In particular, the mapping J~ : O i-» is called Choi- Jamiolkow ski isomorphism. 

Proposition 3.2. For a quantum channel <!>, its dynamical matrix enjoy the properties that follow: 
(i) Dj, = Do,; 
(hJ D<d > 0; 

(mj fr A (D ) - /a, 7HDa>) = 

(iv) |Lo>» - (/®5 ®/)|D a> »;«L (1 ,|L ¥ » - ((D^D^)); (Lo, L«?> - (D^Dy); 

(v) <<B(X), F> - (Do, F 8 X*>/or any X, Y. 

Proof. Write D<d = D = [D,w\. 

(i) D+ = [ZJgjt = [DM* = [D*„ v ] - [Z£J - [L™1 - [DjjJ - D. 



(ii) Let \z) = l£ v=1 z nv \nv)). Then <z| - £^ =1 z^«m/i|. Hence 



flV 



(Z\D\Z) - Xl„ = i Z*D, z Znv 



- 2im=i IW") i s cai l e( l a maximally entangled state. So we have 
|/»«/| = Z^i W){w\ = Z^ v=1 ImXvI ® Im)<v|. 
Since O is completely positive map, ®®Idyt > 0(Vnon-negative integer k), in particular, (0®Id#)(|/)X(/|) > 
0, we get that 

o < (z|(*®idy)G/»«/l)|z> 

N 

= Z <z|[<W<v|) ® \p)(v\]\z) 

fl,V=l 

N N 

= Z Z 4a^<™l$(lj«Xv|)|n> • <a|ju><yl8> 

jU,v=l m,a,njS=\ 



/j Zm/iZnvL™ - Z Z *mfi D Z Znv - a 



flV 

m,^i,n, v=I m,/j,«,v=l 

<I> d>®Id/v 

Obviously, there is an identity in the proof: if &CH B ) — > @CHa), then @CH B ® *Hb) — > @CH A ® 
(3.22) Do = (O ® Id JV )(|/»«/|), 0(p) = TrxtDoC/A ® p l )l 
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Notes: If X = \p)(v\, then |X» - |jU>|v> = \pv), from which it follows that |<B(X)» = L\X)) = L|juv> - 
ZjVi l''7><01^> - Zf,-i £y IU>- Therefore, 0(X) - Ef. 1 I0O1, and 

N 

(m|0)(^)<v|)|n) - <m|0(X)|n>= Y L tf <m|i>0"|n> 

AT 

Since p^ = (p\p\v) = tr(p\v)(p\) = trfl/iXvlp 1 ), we have: 

0(p) = ^ p^vO(|Ai><v|) = J] a>(lpXv|)tr(|juXv^) - ^ Tr A (0(|^>(v|) ® (IpXvlp 1 )) 

= Tr A [(0 ® Id Af )(|/ B »«/ B |)(/ A ® p 1 )] - Tr A [D(p(7 A ® p 1 )]. 

(hi) Since D = [D^], where = [D**], tr A D = [trD^]. Because trD^ v - H^ =1 D, = X^ =1 L™, = 

5^, thus we have tr A D = [6^] = 7 B . Furthermore, Tr(D) = N is trivially. 

(iv) By the operator-sum representation theorem, we have 0(p) = HyTypTt, thus 

l^» = 1 2 O ® r}» = ^ |r ; - ® r*» = ^(/ ® 5 ® /)|r y »|r}» 

j j j 

= ® 5 ® 7)||r y »«r y |» = (/ ® 5 ® /)| 2 ^ ® r})» 

7 y 

= (7 ® S ® 7)|^(2 r y ® r*)» - (7 ® 5 ® 7)|ft(L)» 

i 

= (7®5®/)|D ». 

Therefore «La.|L>p» - «Do>|(7 ® 5 ® 7) 2 |Dvy» - ((D^Dy)); that is (L$, L ¥ > = <Do>, 7>p>- 
(v) 

<y,o(X)> = «y|0(X)» - «r|L*|x» - Tr[L <D |x»«y|] 

= Tr[(|y»«X|) t L (1) ] = <|F»«X|, Lo) = <^(|F»«X|), K(L*)) 
= (Y®X*,D®). 

□ 

For any quantum channel CD, it induces its dual channel & in the following sense: 

(<l>(p), <x) = (p, O^(cr)) for any states p and o~. 
If a CP map is given by the Kraus form <I>(p) = Y*j ^jpfj, then the dual maps reads O^(cr) = £ j ^jpfj. 
Therefore, we have the following proposition into which the most useful results are summarized: 

Proposition 3.3. (i) L* = Y,j Fy ® T*, or D® = £ y |r i »«r y | for <&(•) = ^ Ty • rj. 
( ii) If <t> + is ?/ze <iwa/ channel of a quantum channel O, L^t -To T(L^) - L^, or D^t = 

( Hi) LrQ+sy = rL<$> + sLy, or D r ® +S y - rD® + sDy. 

( iv) the composition Oo^ of two maps O and *P means that L® y = Lq,Ly, or D^oT = ^?(^?(7)$)^?(£)>j/)). 
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(v) L To ® - T r (L®), orD To ® - T A {D®); L<$, oT - T C {L®), orD$, oT - T B (7)<i>). 

(vi) Lrovor = T(U) = L%, orD To(boT = T(D ) = D% = D%. 

Proof, (i) Since L |Z» = |0(X)» = \ TjXT])) = O ® = Tj ® T*. D$ - ft(Lo) - 

Z^(r,®rp = i:,.|r,»«r,|. 

(ii) Obviously, = 2,Tj ® I* - Q^T; ® T*) 1 " - Lj,. Thus it follows from (3) of Proposition 2.9. 
that 

- K(t 4 ,)^(4)^of(L;) 

= [ft o «r(L*)r - cf ° - [W©)]' - Wo) = ^ T(Ad). 

(iii) It is trivially because 

W.vt|X» - KrO + = r|(D(X)» + 

= rLo>|X» + sL^X)) - (rL + sL T )|X». 

D r ®+sW - + ^Dvj/ holds since the reshuffle transformation is linear. 

(iv) L® oW \X)) = |0 o Y(X)» = L^X))) - L^LvplX)). This implies that = K(K(D<d)K(Dy)). 

(v) - 7,t~L<j> = SL^ = T r (L®); similarly, L<j> r = L^Lj- - L®S - T C (L$>). Thus 
D ro<t - ^(Lroo) = ft o 7>(L<d) = KoT r o <R{D®) - 7"a(Ad) 

and 

D<i, or - ft(Lfl>o7-) = ft F c (L<i,) = KoT c o ft(Do,) - 7" b (D<d). 

(vi) L To <$> oT - L T L<$>L T - SL®S - T(L<s,). Thus 

Dro^or = ftCLro^or) - ft o !F(L<d) = ft ° r ° ft(D©) - 7"(D<d) = D$> = D^. □ 

Proposition 3.4. For fwo quantum operations O, *F ora the N -dimensional identical subsystems, "Ha^b 
of a bipartite quantum system 'Ha ® respectively. Then: 

(3.23) = (7 ® S ® 7)(L© ® L T )(7 ® S ® 7). 

TVoo/ p = [p^y] = 1!^v = iP/l*v ® Ia*)(v|, where p^ v - [p»v\, is a N x N block density matrix whose 
entries being N xN scalar matrices. Since 

AT 

(/ ® S ® 7)|p» = (/ ® 5 ® 7) J] |p pv ® |p><v|» 

N 

= J] (7 ® 5 ® 7) 2 |p, v » ® |pv> 

/i,V=l 

/V 

= ^ IP/xv)) ® 
/i,v=l 
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we can get that 

N 

W|P» = KO® l P)(p)»= £ |®(Ppv)®^><v|)» 

//,v=l 

AT 

- ^ (/ ® 5 ® imip^y))) ® ^(I/XXVI)))] 

/U,V=1 

AT 

- J] (/ ® 5 ® 7)[L<D|ppv» ® Lvpl^y)] 

yU,V=l 

AT 

- J] (/ ® 5 ® 7)(Lo> ® ^)[|Ppv» ® lA*v>] 

/U,V=1 

AT 

- (/ ® S ® 7)(Lo> ® Lvp)[ ^ |p pv » ® |/iv>] 

//,v=l 

= (7 ® 5 ® 7)(Lo ® 7>f)(7 ® 5 ® 7)|p». 

□ 

Proposition 3.5. Le? <1>, *P fte ftvo quantum operations on 'Hn- If P,o~ are states in Hn ® < W;v aw <^ 
cr — (O ® *F)(p), 

(3.24) ft(tr) - Un(p)L%. 

Proof. 

|ft(cr)» - (/ ® 5 ® 7)|cr» = (/ ® S ® 7)|(0 ® Y)(p)» = (/ ® 5 ® I)L mw \p)) 
= (I ® 5 ® /)(/ ® 5 ® /XL® ® Ly)(7 ® 5 ® 7)|p» - (L<d ® L ¥ )|^(p)» 
= |Loft(p)4». 

□ 

Lemma 3.6. The composition of two completely positive linear super-operators CD o *F is again com- 
pletely positive. 

Proof. To see that O o *F is completely positive, it suffices to show (O o *P) ® Id^ is positive for any 
iel Obviously, (O o *P) ® \d k - (O ® Id fc ) o (*F ® Id^). Since <D and *F are CP maps, O ® Id^ and 
*F ® Idyt are positive for any JeN, which implies that the composition (O ® Id^) ° OP ® Id/;) is positive 
for any iel □ 

Corollary 3.7. G/vera ftvo Hermitian matrices A,B of the same size N 2 x N 2 . If A,B > 0, then 
R(R(A)R(B)) > 0. 

Proof. We can consider two non-negative matrices A and B of the same size N 2 X N 2 as the dynamical 
matrices for two linear super-operators ®a and <1>b, both acting from M,\< to Mjy, respectively. Now 
A, B > imply that <J>a and ®# are CP maps. Thus their composition ®a ° ®b is CP map by Lemma 
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3.6. It follows from this that, for <I> A o cp e , its dynamical matrix D<s> A o<t> B = < R( ( R(A) ( R(B)) is non- 
negative. The result that follows immediately. Another complicated proof on the present corollary can 
be found in O. □ 

Corollary 3.8. Given a finite set of Hermitian matrices {Dj : j — 1, . . . , n] of the same size N 2 x N . 
IfDj > Ofor all j, then 'R('R(D n ) < R(D n . l ) ■ ■ ■ ft(DO) > 0. 

Proof. For each positive matrix Dj of size N 2 x N 2 , linear super-operator O y determined by Dj is 
completely positive. Thus the composition of n completely positive linear super-operators {Dj : j = 
1, . . . ,n} is denoted by O = <f> n o • • • o Oj. Therefore the dynamical matrix D® for O is equal to 
( R('R(D n )'R(D n ^i) ■ ■ ■ < R{D\)). Since composition preserves completely positivity by the above lemma, 
cD is completely positive, therefore > 0. □ 

Proposition 3.9. The Hilbert-Schmidt inner product, i.e., (X, Y) = Tr{X^Y), on the matrix space Mn 
induces another inner product in the space of linear maps ££(Mn, Mn)- 

Proof. Let [E a : a = I,... ,N 2 } and {F a : a - 1, . . . ,N 2 } be orthonormal bases in Mn, where 
(E a , Ep) = (F a , Fp) = 8 a p. We need only to prove that 

YS=x Tr<S>{E a )^{E a ) = TrO(F ff )t«F (Fa ). 
Since \E a )) = j£ x \Fp)){{Fp\E a )) = Ylpt x {{Fp\E a ))\Fp)), E a = ^((Fptf^Fp. 

N 2 N 2 N 2 

2 TrtKE a ?V(E a ) = Z Z «^8l^»«^l^»TrO(F i8 ) t T(F y ) 

a=\ a= \ P,y=l 

N 2 N 2 

= Z Z «^rl^»«^l^»TrO(F j8 ) t ^(F r ) 

a=l p,y=l 
N 2 N 2 

= Z « F rl(Z l^»«^l)l^e»Tra>(F^) t Y(F r ) 

p t y=\ a =\ 
N 2 N 2 

= {{F y \Fp))TrQ>{Fp)^{F y ) = £ 5 /Jy TrcD(F i8 ) t ^(F r ) 

p,y=l p,y=\ 
N 2 

= ^ Tr<D(F a ) tl F(F ff ). 



a=\ 



□ 



Now we define the inner product of two linear super-operators O and *F (see [2]) as follows: 



(3.25) (O, «P) = J] imE a )^{E a ). 

a=\ 

Using this correspondence it is possible to introduce two different bases, associated to the bases 
[E a } N =l ,{Fp}^ v in the space of linear maps: 
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1) Type-I basis {A^} in S£(Mn, Mn) is defined by 

(3.26) A af3 (X) = E a (Ffj, X) = E a TrF^X, X e M N ; 
and 

2) Type-II basis {& a p} in ££ (Mn, Mn) is defined by 

(3.27) ® af} (X) = E a XFp, X € M N . 

Indeed, 1) Let Z^g = i c a p& a p = for some scalars c a p € C. This implies that E^s=i CapA a p(X) = 0, 
in particular, for X = F y (y = 1, . . . , A^ 2 ), we have: 

Since {£"0,} is linearly independent, = 0(a, y = l,...,N 2 ). We have also that (A a p,A MV ) = 
2S,/=i Tr[A ff/? (|/>01)^ v (|/c></|)] - ^fyv Furthermore, L Aff/J = |£«»«F / ,|. 

2) Let 11^8=1 c a /j® a /j = for some scalars c a p e C. This implies that 2^8=1 c ap®ap(X) = 0> we have: 



- «=i C aj 8® aj g(X) - 1<a£=l c apE a XFp, 



which means that 



N 2 N 2 N 2 

afi=\ afi=\ afi=l 



^ c al3 E a ®F*p 



\X)); 



afi=\ 



i.e., ^j=i CapE a ®F*p = since X is arbitrary. Because of the independence of the set {E a ® F^„ v 
this implies that c a p = 0(a,fi = I,... ,N 2 ). And we have also that (& a p, ©^y) = 5 ail 8p v . Furthermore, 
L &ap = E a ® F*.D 

Remark 3.10. Therefore, according to two kind of the above-mentioned bases, we can expanding 
any mapping <1> € j£f (Mn, Mn) with respect to Type-I and Type-II bases, respectively, to get two 
expressions that follow: 

N 2 N 2 

(3.28) 0> = PapKp = ^ q a /3®a/3- 

a,p=l aJ3=l 

Now Lo = xTj3=i Pop\E a ))((Fp\ = Sj^i q ap E a ® F*. We write P = [p afj ], Q = [q ap \. 

There is natural question to be asked: what is the relationships among these matrices P, Ql (see iflQl ) 

Proposition 3.11. With the above notations, 

(3.29) (K p , % v ) = (® ap , A MV ) = TKElE^Fl). 
Thus 

(i) Pap = EjJ=i Tr(EiE^FpFl)q^; 
(U) q a p - Z^ v= i TriElEftFfiF^Pfiy. 
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Proof. By the definition of the inner product in the space of linear maps, 

N N 

(Kp,®^) = J] TrCCA^dOOD^Q^vdOOl)) = j <i|F^)Tr(4^|j)<7l4) 

ij=i t,y=i 

AT Af 

= £ < W> • </Fj£^|0 = J^(i\F fi FtE%\i) = iTiFpFlElE^) 

1,7=1 i=l 

- TviElE^FpFl). 
Similarly, we have also: (& a p, A^y) = Tr(E a E M FpFl). Since 



AT 



(U,V=1 



AT- 



2 



jU,v=l 



1) and 2) is trivial. 



Remark 3.12. (i) A special case is provided by the choice E a = F a (or E a - F a - where 
are an orthonormal basis for C ) (see |[TTl). 

(ii) Since |/» = \Zi\m) = Zil"), |/»«/| - Xijmm = Z y l0Ol®l*X/l- We know that 
I® I - Zff=i |Fa)X(F<*l when {£' Q .}^ 1 is an orthonormal basis for Mn- Thus we have: |/))((/| = 
<R{I ® I) = Z^Ii K(\E a ))((E a \) = Z^Ii E a ® E* . If there is another orthonormal basis {Fp}^ 
for M N , we still have: |/»«/| = F fi ® F^. Therefore, Z^Ij E a ®E* a = Z^ F fi ®F* = 
Zf 7= i I0O1 ® I 1 ') 01 = Furthermore, we have the swap operator 5 - Zf/=i 10)0*1 - 
l£ tl E a 9Et = X* l F fi *Fl. 

(iii) In fact, given two orthonormal bases [E a }^^ and [F a }^_ 1 in AIa?, the relation 

N 2 

(3.30) if (Mat, M w ) 3 <D — > A© = J] 0(F ff ) 8 F ff e M w ® M w 

a=i 

defines an isomorphism between J£{Mn, Mn) and Mn® Mn- The isomorphism is an isometry: 

N 2 N 2 N 2 

(Ac, A>p> = <J] ®(E a ) ® F a , J] ® Fg) = ^ ® F a , V(E fi ) ® Fp) 

a=\ jS=l afj8=l 

W 2 A? 2 

afi=\ a, (3=1 
N 2 

= ^<0)(F a ), y(E a )) = <<D, *F>; i.e., (A^, A>p> - <<D, ¥>. 

a=l 
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4. Best separable approximation for states 

In this section we recall the so-called optimal and the best separability approximation(OSA and 
BSA respectively). Although the results below have been proven in (HlH], we give the framework for 
our convenience. Other results involved can be found in [14j. In the method of BSA, for any density 
matrix p there exist a "optimal" separable matrix p* and "optimal" non-negative scalar A such that 
p - Ap* > 0. We describe these results involved that follow: 

Definition 4.1. A non-negative parameter A is called maximal with respect to a (not necessarily nor- 
malized) density matrix p, and the projection operator P = \ifr)(iff\ if p - AP > 0, and for every e > 0, 
the matrix p - (A + e)P is not positive definite. 

Definition 4.2. A pair of non-negative (Ai,A2) is called maximal with respect to p and a pair of 
projection operators Pi = \if/i)(if/\\, P2 = I1A2X1A2I, ifp — AiPj - A2P2 ^ 0, Ai is maximal with respect 
to p - A2P2 and to the projector P\,A2 is maximal with respect to p - A1P1 and to the projector P2 
and and the sum Ai + A2 is maximal. 

Theorem 4.3. For any density matrix p (separable, or not) and for any (fixed) countable set V of 
product vectors belonging to the range of p, there exist A(V) > and a separable matrix 

a 

where each projector P a is generated by some product vector in R(p), and all A a > 0, such that 
8p = p - Ap* > 0, and that p* s (V) provides the optimal separable approximation (OSA) to p since 
Tr(Sp) is minimal or, equivalently, A is maximal. There exists also the best separable approximation 
p* for which A - maxy A(V). Obviously, A(V) < A(V') when V c V. 

Theorem 4.4. Given the set V of product vectors in the range R(p) of p, the matrix p* - A a P a is 
the optimal separable approximation OSA) of p if: 

1 ) all A a are maximal with respect to p a — p — Y^a+a A- a <P a ', and to the projector P a >; 

2) all pairs (A a ,Ap) are maximal with respect to p a p — p — Y^a+afi A ff ' F ' a i , and to the projection 
operators (P a , Pp). 

Theorem 4.5. (The uniqueness of the BSA ) Any density matrix p has a unique decomposition p = 
Ap s + (1 - A)5p, where p s is a separable density matrix, 6p is a inseparable matrix with no product 
vectors in its range, and A is maximal. 

5. Best separable approximation for operations 

We cab define separable CPM; that is, <1> is separable if its action can be expressed in the form 

®(p) = YU{Sk®T k )p(Si®T k )\ 
for some integer n and where St and 7^ are operators acting on 'Ha/b, respectively. Otherwise, we 
say that it is nonseparable. Up to proportionality constant, separable maps are those that can be im- 
plemented using local operations and classical communication only. 
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Let us consider two systems, A and B, spatially separated, each of them composed of two particles 
(A12, and 61,2)- Let us consider a CPM O acting on systems A\ and B\, where <&(p) = Y*k M (i) p(M w ) t 
and M (k) acting on < H A , ® 1-(b { , where M ik) = S k ® Tk for each index k. Now <1) induced another super- 
operator acting on 1-l A ® *Ke in the following sense: 

= Z*CS* ® f*)X(Sl® f t )t, 
where 5 fc = ® id (A2) and f k = T k ® id (B2 \ and X acting on 74 ® 

We are interested in whether this CPM can create "nonlocal" entanglement between the systems A and 
B. We define the operator E Ai a 2 ,BiB 2 acting on ( H A ® ( H B , where 1-t A - < Ha 1 ® < Ha 2 an d 'Kb = 'Kbi ®*Hb 2 , 
and dimCHii) - dimCWfl.) - J, as follows: 

£A a A 2>Bl B 2 - (® (AlBl) ®id (A2B2) )(PA I A 2 ®Pl» 1 B 2 ) = ^(PA 1 A 2 9PB 1 B 2 ) = Z*(5*®^t)(i , A,A 2 ®PB I B 2 )(iS*® 
^t) 1 " = Y> k $kPA l A 2 T k h®{T k PB l B 2 Tk). 

Here, /% A2 - |f> AlA2 (f | with \W) Ai a 2 = ^= 2^ =1 |m) A] ® |m> Az , and P BlB2 = |^> Sl s 2 <^| with 
l vp >BiB 2 = -^f I^!=i ly«>Bi ® ly">B 2 ' where d m > : m = l,---,d) and {|//> : ju = 1, . . . , d} are an or- 
thonormal basis for < H Ay j Al and < Hb 1 /b 2 respectively. The map <I> is understood to act as the identity on 
the operators acting on "H A2 and 'Hb^ The operator E has a clear interpretation since it is proportional 
to the density operator resulting from the operation <1> on systems A\ and B\ when both of them are 
prepared in a maximally entangled state with two ancillary systems, respectively. E is called Choi ma- 
trix for the bipartite super-operator <1>, or the mapping O — > £"(<!>) is called Jamiolkowski isomorphism 
for the bipartite super-operator O. 

Now in general for 0(p) = ^M^pCM®) 1 ", where M« = Y>mn,nv M mlnv\m){n\ ® |/x)(y|, then <D 
induced another super-operator like above as follows: 

5>(X) - Y,k^ k)x (^) , 
where = 2 mn>/iV M® ;/ , v |^Xn| ® and |^>^| = |m)(n| ® id (Az) and |^><^| - |ju>(v| ® id (B2) . 

E AlAM = (3> (AlBl) ® id (A2fi2) )(P AlA2 ® P BlB2 ) 

= $(P AlA2 ® P Bl B 2 ) = J] M^\P Al A 2 ® P Bl B 2 ) (m^) + 

/ i / 1 lyl mn,M lyl m'n',ti'v' s 
k mm'nn' ,/i/i'vv' 

(\m)(n\ ® id (A2) ) J P AlA2 (|m'>(«'| ® id^) 1 " ® (|ju><v| ® id (B2) )P BlB2 (|//><v'| ® id <2fe) ) f 
-V V M W \M (k) V 

k mm! nn' ,im' vv' 

(\m)(n\ ® id (A2) ® |/x)(v| ® id (B2) )(P A]A2 ® / J Bl B 2 )(|m , )<n'| ® id (Az) ® \n'){v'\ ® id (B2) ) f . 

If we define vec(M w ) = Z m „ >iUV M { ^ v \m)\n)®\n)\v) - Y.mn^v M^ v \mn)®\nv) - Y.mn^y M^jJmnfiv), 
then Zk vec(M«)vec(M«)t = £ AiA2 , b ,b 2 . If M« = A k ® then vec(M«) = |A*»|£*». 

Proposition 5.1. If<t> is a quantum operation on a bipartite quantum system, then O is separable if 
and only if its dynamical matrix D<d is separable. 
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Proof. By the definition of separable quantum operation, cp(p) = £ ; (Aj ® B,)p(A, ® B;)' when O is 
separable. Now the dynamical matrix for the separable operation <& is = vec(A,- ® fi,)vec(A, ® 

Definition 5.2. Given quantum operation <I> on bipartite quantum system 'Hi ® "7^2 with dimTYi = 
dimT-^ = A 7 ", the dynamical matrix D$ for cD can decomposed as = AD S + (1 - A)D e in terms of 
the BSA decomposition for state. Then the separable operation $>bsa determined by AD S is called best 
separable operation approximation for O. Q>ent = ® — ®bsa is called pwre entanglement-produced 
operation part for O. 

Remark 5.3. If there is another decomposition = D' s + D' e for which D' s is just separable, then: 
AD S - D' s > by the uniqueness of the BSA. Thus the decomposition cD = ®bsa + ®ent is unique. 

By operator-sum representation theorem, <l>(p) = TjiefFjpFj = T^jeG^ jpGj, where max(|F|,|G|) < 
N 4 . Let 

I = {i e F : F { = A t ® 5,}, J - { j e G : Gj = C,- ® D ; }. 
Write Y(p) - Zm^J and Y(p) = ^GjpG^Tfp) = Z ieF \l and ¥'(p) - Z yeOyrG/pGt. 

Theorem 5.4. T = ¥ = O bsa . 

Proof. Apparently, D® = £)y + Df> , where Dx is separable since T is separable operation. Hence 
it follows from the uniqueness of the BSA that AD S - Dx > 0. If, otherwise, AD S - Dx > 0, then 
d ®bsa-v ~ D<s>BSA~ D r > 0, that is, ® BSA -Y is CP map and separable, soD$ = [Dr+Do BSA -r]+[Dr'- 
D,d bsa -y], where Dx' - Dcd bsa _y > 0, contradict with the fact that there is no factorizing operational 
element for Y'. Therefore AD S - Dx = 0, equivalently, Y = Qbsa- The theorem is proved. □ 
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